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The paper analyzes the phenomenon ol turbulent thermal diffusion in the Earth atmosphere, its 
relation to the turbulent diffusion and its potential impact on aerosol distribution. This phenomenon 
was predicted theoretically more than 10 years ago and detected recently in the laboratory experi- 
ments. This effect causes a non-diffusive flux of aerosols in the direction of the heat flux and results 
in formation of long-living aerosol layers in the vicinity of temperature inversions. We applied the 
theory of turbulent thermal diffusion to the GOMOS aerosol observations near the tropopause in 
order to explain the shape of aerosol vertical profiles with elevated concentrations located almost 
symmetrically with respect to temperature profile. We demonstrate that this theory is in good 
agreement with the observed profiles of aerosol concentration and temperature in the vicinity of 
the tropopause. In combination with the derived expression for the dependence of the turbulent 
thermal diffusion ratio on the turbulent diffusion, these measurements yield an independent method 
for determining the coefficient of turbulent diffusion at the tropopause. We also derived a practi- 
cally applicable formulation for dispersion of atmospheric trace species which takes into account the 
phenomenon of turbulent thermal diffusion. We evaluated the impact of turbulent thermal diffusion 
to the lower-troposphere vertical profiles of aerosol concentration by means of numerical dispersion 
modelling, and found a regular upward forcing of aerosols with coarse particles affected stronger 
than fine aerosols. 

PACS numbers: 



I. INTRODUCTION 

Various aspects of turbulent diffusion of aerosols in 
the atmospheric flows have been extensively investigated 
in the past [e.g., Monin and Yaglom, 1975; Csanady, 
1980; Maxey, 1987; Flagan and Seinfeld, 1988; Wyn- 
gaard, 1992; Fessler et al., 1994; Blackadar, 1997]. In 
particular, the turbulent diffusion (eddy diffusivity) has 
been comprehensively studied for the low-order closures 
of the turbulent dispersion equation. However, certain 
important features of turbulent transport of aerosols in 
stratified flows have been found only recently. In partic- 
ular, a new phenomenon of turbulent thermal diffusion 
has been predicted theoretically by Elperin et al. [1996, 
1997a] and detected in the laboratory experiments in sta- 
bly and unstably stratified turbulent flows by Eidelman 
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et al. [2004, 2006a, 2006b] and Buchholtz et al. [2004]. 
The phenomenon of turbulent thermal diffusion (TTD) in 
turbulent stratified flows results in the non-diffusive flux 
of aerosols and gaseous admixtures in the direction of the 
heat flux. Particles are accumulated in the vicinity of the 
minimum of the mean temperature of the surrounding 
fluid. This phenomenon causes formation of large-scale 
inhomogeneities in spatial distribution of aerosol parti- 
cles in the vicinity of temperature inversions. 

The effect of turbulent thermal diffusion has been de- 
tected in two experimental set-ups: oscillating-grids tur- 
bulence generator [Eidelman et al., 2004, 2006a; Buch- 
holtz et al., 2004] and multi-fan turbulence generator 
[Eidelman ct al. , 2006b] . The experiments have been per- 
formed for stably and unstably stratified fluid. In these 
experiments, even with strongly inhomogeneous temper- 
ature fields, particles in turbulent fluid accumulate in the 
regions of temperature minima, in a very good agreement 
with the theory of turbulent thermal diffusion. 

In spite of the previous comprehensive theoretical and 
laboratory studies of the phenomenon of turbulent ther- 
mal diffusion, the observational evidence and quantita- 
tive evaluation of its importance in the Earth atmosphere 
have not been investigated until now. In the present 
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study, we analyzed the GOMOS observations in the vicin- 
ity of the tropopause and explained the shape of aerosol 
vertical profiles with elevated concentrations near the 
minimum of temperature. The contribution of the effect 
of turbulent thermal diffusion to the lower-troposphere 
vertical profiles of the aerosol concentration was investi- 
gated via aerosol dispersion modelling. 

The existing theory of turbulent thermal diffusion 
[Elperin et al., 1996, 1997a, 1997b, 1998, 2000a, 2000b, 
2000c, 2001] does not take into account the structure of 
the atmospheric stratified flows. In this paper, the theo- 
retical approach of Zilitinkevich et al. [2007, 2008] is fur- 
ther developed and applied to study the effect of stratifi- 
cation on turbulent transport of aerosols. To this end we 
derive the budget equation for the turbulent flux of parti- 
cles in stably stratified flow. This allows determining the 
dependence of the turbulent diffusion and turbulent ther- 
mal diffusion coefficients on the flux Richardson number. 
We demonstrate that the coefficients of turbulent ther- 
mal diffusion and turbulent diffusion decrease with the 
increase of the flux Richardson number. 



II. TURBULENT THERMAL DIFFUSION AND 
TURBULENT FLUX OF AEROSOLS 

A. Mechanism of turbulent thermal diffusion 

Let us discuss the physics of the phenomenon of tur- 
bulent thermal diffusion. We consider inertial particles 
(aerosols) suspended in the turbulent fluid flow with large 
Reynolds numbers. Particle concentration n p = N + n 
is characterized by the mean value, N, and fluctuations, 
n (measured in m -3 ). Evolution of the number density 
n p (t, r) of small inertial particles in a turbulent flow is 
determined by the following equation: 



dn p 



V-(n p v) =DV 2 



(1) 



where v is a random velocity field of the particles which 
they acquire in a turbulent fluid velocity field u, and D 
is the coefficient of molecular (Brownian) diffusion. We 
assume here for simplicity that the mean velocity is zero, 
and we do not take into account the effect of particles 
upon the carrying fluid flow. The velocity of particles 
v depends on the velocity of the surrounding fluid and 
it can be determined from the equation of motion for a 
particle. When p p 3> p, this equation represents a bal- 
ance of particle inertia with the fluid drag force produced 
by the motion of the particle relative to the surrounding 
fluid, dv/dt = (u — v)/t s , where t s is the particle Stokes 
time, p is the fluid density and p p is the material density 
of a particle. Solution of the equation of motion for small 
particles yields: 



[see, e.g., Maxey, 1987]. The second term in Eq. ([2]) de- 
scribes the difference between the local fluid velocity and 
particle velocity arising due to the small but finite iner- 
tia of the particle. In this study we consider low Mach 
numbers turbulent flow with V- u = — p^ 1 (u-V)p 7^ 0. 
Equation |2]) for the velocity of particles and Navier- 
Stokes equation for the fluid for large Reynolds numbers 
yield 



V-v 



Vu 



1 
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= --( u -V)p+^V^ + 0(r s 2 ) 
P P 



(3) 



V = U - T, 



- + (u.V)u 



+ 0(r s 2 ) , 



(2) 



where p is the fluid pressure. 

The physical mechanism of the phenomenon of tur- 
bulent thermal diffusion for inertial particles can be ex- 
plained as follows. Due to inertia, particles inside the tur- 
bulent eddies drift out to the boundary regions between 
the eddies (the regions with the decreased velocity of the 
turbulent fluid flow). Neglecting non-stationarity and 
molecular viscosity, the estimate based on the Bernoulli's 
law implies that these are the regions with the increased 
pressure of the surrounding fluid. Consequently parti- 
cles are accumulated in the regions with the maximum 
pressure of the turbulent fluid. Indeed, due to the iner- 
tia effect V- v oc (r s /p) V 2 p =/= even for incompressible 
fluid flow [see Eq. ((3)]. On the other hand, for large 
Peclet numbers, when we can neglect the molecular dif- 
fusion of particles in Eq. ([1]), V-v oc —dn p /dt. This 
implies that in regions with maximum pressure of tur- 
bulent fluid (i.e., where V 2 p < 0) there is accumulation 
of inertial particles (i.e., dn p /dt cx — (r s / p) V 2 p > 0). 
Similarly, there is an outflow of inertial particles from re- 
gions with minimum pressure of fluid. Note that in cloud 
physics, the effect of local accumulation of particles be- 
tween turbulent eddies (in the regions with maximum 
fluid pressure) has been invoked extensively in order to 
elucidate the mechanism of rain formation [see, e.g., Pin- 
sky and Khain, 1997; Shaw, 2003; Collins and Keswani, 
2004; Wang et al., 2006; Khain et al., 2007]. 

In case of homogeneous and isotropic turbulence with- 
out external large-scale gradients of temperature, a drift 
from regions with increased (decreased) concentration of 
particles by a turbulent flow of fluid is equiprobable in 
all directions. Therefore pressure (temperature) of the 
fluid is not correlated with the turbulent velocity field 
and there exists only turbulent diffusion of particles. 

Situation drastically changes in a turbulent fluid with 
a mean temperature gradient. In this case, the heat flux 
(u9) is not zero, i.e., fluctuations of fluid temperature, 
8, and velocity are correlated. We consider low-Mach- 
number flows (M = u/c s -C 1, c s is the sound speed) and 
study mean-field effects. For low-Mach-number flows, the 
mean fluid mass flux (u p') is very small (~ 0(M 2 )) [see, 
e.g., Chassaing et al., 2002], i.e., the fluctuations of the 
fluid density p' and velocity u are weakly correlated. On 
the other hand, fluctuations of pressure must be corre- 
lated with the fluctuations of velocity due to a non-zero 
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turbulent heat flux, (u9) ^ 0. Indeed, using the equation 
of state for an ideal gas we find that 



P_ 

P 



EL 
p 



T ' 



(4) 



(see below) and (up}/P = (u0)/T, where P, T and p 
are the mean fluid pressure, temperature and density, 
respectively. Therefore, the fluctuations of temperature 
and pressure are correlated and the pressure fluctuations 
cause fluctuations of the number density of particles. In- 
deed, increase (decrease) of the pressure of surrounding 
fluid is accompanied by accumulation (outflow) of the 
particles, respectively. The direction of the mean flux 
of particles coincides with the direction of the heat flux 
of temperature - towards the minimum of the mean tem- 
perature. Therefore, the particles are accumulated in this 
region [for more details, see Elperin et al., 1996, 1997a]. 

Equation (j4]) is obtained as follows. Averaging the 
ideal gas equation yields equation for the mean pressure: 
P = (Kfj/m^) (pT + (p' 9)), whereas the equation for the 
pressure fluctuation is p = (ftfe/m^) (p9 + p' T + p' 9 — 
(p' 6}). Here Kb is the Boltzmann constant and is the 
mass of molecules of surrounding fluid. Nonlinear terms 
p' 9 and (p' 9) in this equations are of the order of 0(M 4 ) 
[see, e.g., Chassaing et al., 2002], and for fluid flows with 
small Mach numbers, they can be neglected. This yields 
Eq. (g]) for the ratio p/P. 

In order to demonstrate that the directions of the mean 
flux of particles and the turbulent heat flux of tempera- 
ture coincide, let us consider fluid with the mean temper- 
ature gradient (Fig. 1). Assume that the mean temper- 
ature T 2 at point 2 is larger than the mean temperature 
Ti at point 1. Let us consider two small volumes a and 
b located between these two points and let the direction 
of the local turbulent velocity of the volume a at some 
instant be the same as the direction of the turbulent heat 
flux, (u9), i.e., towards the point 1. Let the local turbu- 
lent velocity of the volume b be directed at this instant 
opposite to the turbulent heat flux (i.e., to the point 2). 
In a fluid with an imposed mean temperature gradient, 
fluctuations of temperature 9 and velocity u are corre- 
lated. Positive temperature fluctuations result in positive 
pressure fluctuations. Consequently, the fluctuations of 
the temperature 9 and pressure p are positive inside the 
volume a and negative inside the volume b. The fluctua- 
tions of the particle number density n are positive in the 
control volume a (because particles are locally accumu- 
lated in the vicinity of the maximum of pressure), and 
they are negative at the volume b (because there is an 
outflow of particles from regions with a low pressure). 
Consequently, the mean flux of particles is positive in 
the volume a (i.e., it is directed to the point 1), and it 
is also positive inside the volume b (because fluctuations 
of velocity and number density of particles are negative 
in the volume b). Therefore, the mean flux of particles is 
non-zero and directed, as well as the turbulent heat flux 
(u9) of temperature, towards the point 1. 
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FIG. 1: Mechanism of turbulent thermal diffusion. 

B. Turbulent flux of aerosols in stratified fluid 

The theory of turbulent thermal diffusion developed 
previously [see Elperin et al., 1996, 1997a, 1998, 2000c, 
2001], does not take into account the effects of the 
Richardson number and anisotropy of turbulence. In this 
study, the parameters of turbulent thermal diffusion are 
derived as functions of the flux Richardson number and 
other turbulence characteristics for stably stratified tur- 
bulent flows. Equation for the evolution of the mean 
number density N of particles reads: 



dN 
~dt 



V- [JV(U + W 9 ) +F (ll) ] = 



(5) 



where U = (U\,U2,Uz) is the mean fluid velocity (e.g., 
the wind velocity) , W g = t s g is the terminal fall velocity 
of particles, g is the acceleration of gravity. The equation 
for the turbulent flux of particles, = (vn), is derived 
in Appendixes A, B, C using different approaches. The 
vertical component of the particle turbulent flux Fj in- 
cludes contributions of turbulent diffusion and turbulent 
thermal diffusion: 



= Vf N -K D V Z A , 



(0) 



where Kp is the coefficient of turbulent diffusion and 
V 2 ee d/dz. The effective velocity V cff caused by TTD 
is given by the following equation: 



V c 



-t T (v V-v) 



(7) 



where ty = i/EJ is the turbulent time, I is the tur- 
bulent length scale, Ek is the turbulent kinetic energy. 
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Equation (Jjj for the effective velocity has been derived 
using different rigorous methods by Elperin et al. [f996, 
f997a, 1998, 2000c, 2001], Pandya and Mashayek [2002] 
and Recks [2005]. Note that even a simple dimensional 
analysis yields the estimate for the effective velocity V 
that coincides with Eq. ((7]). Indeed, let us average Eq. (fT]) 
over the ensemble of the turbulent velocity field and sub- 
tract the obtained averaged equation from Eq. JJJ). This 
yields equation for the fluctuations n of particle number 
density 



rim 

^- + V-(Av + Q) 



DAn 



(8) 



where Q = vn — (vn) is the nonlinear term. The mag- 
nitude of dn/dt + V- Q DV 2 n can be estimated as 
n/tT- Therefore, the turbulent component n of particle 
number density is of the order of n sa — tT V- (Av) = 
-t T [NV- v + (v-V)iV]. Now let us calculate the turbu- 



lent flux of particles F t 



(n) 



(vi n): 



F, 



in) 



-Nt T {Vi V- v) - t T {ViV^VjN , 



(9) 



where the first term in the right hand side of Eq. ((9]) de- 
termines the turbulent flux of particles caused by turbu- 
lent thermal diffusion: —N tr (t>i V • v) = Vf s N, while 
the second term in the right hand side of Eq. (J9j) deter- 
mines the turbulent flux of particles caused by turbu- 
lent diffusion: tx (viVj)W jN = KpViN. In the latter 
estimate we neglected the anisotropy of turbulence for 
simplicity. 

For non-inertial particles advected by a turbulent fluid 
flow, particle velocity v coincides with fluid velocity u, 
and Vu « -(u-V)p/p « (u ■ V)T/T, where p and T 
are the density and temperature of the fluid. This for- 
mula takes into account the equation of state for an ideal 
gas but neglects small gradients of the mean fluid pres- 
sure, i.e., (V zj o)/p w -(V Z T)/T. Thus, the effective ve- 
locity ([7]) for non-inertial particles is determined by the 
following equation: 



v: 



off 



-K D 



T 



(10) 



Alternative derivations of Eqs. © and (fTTJ]) for non- 
inertial particles are also presented in Appendixes A and 
B. 



C. Effective velocity of aerosols caused by 
turbulent thermal diffusion 

Let us now consider incrtial particles and determine 
the dependence of the vertical component of the effective 
velocity V z s on parameters characterizing turbulence, 
the mean temperature profiles and particles. Note that 
the deviation of the particle velocity v from the fluid ve- 
locity u is small, but the deviation of V- v from V- u is 
not small (see Eq. Equations © and J7J) yield 



V 



eff 



-K 



D ■ 



(ii) 



where the second term in Eqs. (jlip describes the contri- 
bution of the particle inertia effect to the effective veloc- 
ity. In order to determine this contribution we use the 
following formulae: 



— (u z V p = — y — 
P T 



(12) 



(u z V 2 9) = — ln(Re)5(Re,a*)V z T, (13) 
o it 

[see Elperin et al., 1996, 1998, 2000a], where 9 are fluc- 
tuations of the temperature, T is the mean tempera- 
ture, Lp 1 = \X7 Z P/P\ and P is the mean fluid pressure, 

1 /2 

Re = tEjl /vis the Reynolds number, v is the kinematic 
viscosity, and a* is the particle size (e.g., the particle di- 
ameter). In Eqs. pip and (fT2")l we neglected small effects 
~ 0(jl ). For derivation of Eq. (fT2")) we took into ac- 
count the equation of state, neglected the flux of fluid 
mass (u p 1 ) for the low-Mach-number flows, and used the 
identity t s — pW g Lp/P, where |V Z P| = pg. Derivation 
of Eq. (fHij) is given in Appendix C. 

When the particle size a* < a cr , the function 
B(Re, a») = 1, and for a* > a cr the function i?(Re, a*) 
is given by B(Re, a*) = 1 — 3 ln(a*/a cr )/ ln(Re) [see 
Elperin et al., 2000a], where the critical particle size is 
a cr = l v (p/ pp) 1 / 2 , and l v = £Re~ 3 ^ 4 is the Kolmogorov 
viscous scale of turbulence. The vertical component of 
the effective velocity caused by turbulent thermal diffu- 
sion can be rewritten as follows: 



Vf = 



J^D &TD „ 

P T 



Kr 



T 



(14) 



where a TD = 1 + Ktd/Kjj. In order to determine the 
explicit form of the coefficient of turbulent diffusion Kd 
and the coefficient Ktd in Eq. (jT4j) , we have to use some 
turbulent closure model. In this study we use the tur- 
bulent closure model for stably stratified flows developed 
by Zilitinkevich et al. [2007]. This model is based on 
the budget equations for the key second moments: tur- 
bulent kinetic and potential energies, and vertical turbu- 
lent fluxes of momentum and buoyancy (proportional to 
potential temperature). This model takes into account 
the non-gradient correction to the traditional buoyancy 
flux formulation and implies the existence of turbulence 
at any gradient Richardson number. Predictions from 
this model are consistent with the available data from 
atmospheric and laboratory experiments, direct numeri- 
cal simulation and large-eddy simulation. 

The turbulent closure model by Zilitinkevich et al. 
[2007] is developed in the geophysical approximation for 
stably stratified flow, whereby the vertical mean fluid ve- 
locity, U3, is negligibly small compared to the horizontal 
velocities, Ui and Ui- A useful approximation for sta- 
bly stratified flows is that the horizontal derivatives of 
the mean velocity components XJ\p are negligibly small 
compared to their vertical derivatives. An effect of the 
horizontal derivatives of the mean velocity components 
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Ui 2 on turbulent transport of particles in stably strati- 
fied flows is a subject of a separate study. 

Using the turbulent closure model by Zilitinkevich et 
al. [2007] we obtain formulas for the coefficient of turbu- 
lent diffusion Kjj: 



K D = 2C n E l J 2 £ z 



1-C 



D 



Ri/ 



2C K A Z (1-Bi f ) 



1 + Cd C n -t— 



(15) 



and for the coefficient Ktd in Eq. (|14p that accounts for 
particle inertia: 



K TD = -C n W g Lp\n(Re)B(Re,a*) 

1 + Cd C n -j— 

E z 



(16) 



Derivation of Eqs. (fT5|) and (jTHJ) is given in Appendixes 
C and D. In the above equations, A z = E z /Ek is 
the vertical anisotropy of turbulence and E z is the ver- 
tical turbulent kinetic energy. The gradient Richard- 
son number is defined as Ri = AT 2 /S 2 , where N 2 — 
(3 d® I dz is the squared Brunt- Valsala frequency, S — 

1 /2 

[{dUx/dz) 2 + (dU 2 /dz) 2 } 1 is the mean velocity shear, 
is the mean potential temperature (or the mean vir- 
tual potential temperature accounting for specific hu- 
midity), that is defined as 6 = T(P*/P) 1_ t _1 . Here 
P* is the reference value of the mean fluid pressure P, 
7 = c p /c v = 1.41 is the specific heat ratio, (3 = g/T* is 
the buoyancy parameter and T* is a reference value of 
the mean temperature. The flux Richardson number is 
defined as Ri/ = —(3F z /KmS 2 , where Fi = (ui9 p ) is the 
flux of the potential temperature, 9 P are the fluctuations 
of the potential temperature, Km is the eddy viscosity, 
E z = E z j (S£ z ) 2 , £ z is the vertical turbulent length scale, 
C ni Ck and Cd are empirical dimensionless coefficients. 

The turbulent closure model suggested by Zilitinkevich 
et al. [2007] yields the following formulas for the param- 
eters in Eqs. ©, (fT4"| - (fll))) : the vertical anisotropy of 
turbulence A z , 



A 7 . 



E, 



2C*-tf r (l-Ri / ) 
the vertical turbulent length scale, 

Ri/ 1 4/3 



1 



Ri^ 



(17) 



(18) 



where * T = 0.2(1 - Ri/), 

2C K ^r 



E z = 



3(1 + C r ) 



[C r * 3 - (a*3 + 3)Ri/] , (19) 



and 5*3 = 1 — C3 Ri/ . For very large gradient Richardson 
numbers, the flux Richardson number Ri/ — » Ri^° = 0.2, 

and the function E z E™ = 2C 9 C K ^Rif, where 



Cg, C3 and C r are empirical dimensionless coefficients. 
Since the coefficient of turbulent diffusion Kd should be 
positive, the empirical constant Cd < 2/3. The empiri- 
cal constants have been determined by comparing results 
from the local closure model by Zilitinkevich et al. [2007] 
with data from laboratory and field experiments, large- 
eddy simulations (LES) and direct numerical simulations 
(DNS): C 3 = 2.3, C r = 3, C K = 1.1, C g = 0.3, C n = 0.3 
and Cd = 0.3. 

The horizontal components of the particle flux are 
given by the following formulas: 



F. 



O) 



"Cr, 



E 



1/2 



-C n 



F: 



(n) 



E 



1/2 ' 



(20) 



where i = 1,2. Equation (|20p implies that in neutral 
stratification the horizontal turbulent flux of particles is 
of the same order as Pj , whereas in strongly stable 
stratification it is approximately by a factor of z/L larger 

(n) 

then F z . This flux generally deviates from the mean 
wind direction, thus contributing to the horizontal cross- 
wind dispersion. 

The steady-state solution of Eq. ([5]) reads: 

V ' A (21) 



N 



T 



"9 







where the turbulent thermal diffusion ratio a TD is given 
by the following equation: 



a, 



Ktd 
~K^ 



1~C D 



W g Lp ln(Re)P(Re,a» 



3P Z 1/2 £ z 



Ri, 



2C K A Z (1-Ri/) 



(22) 



For gases and very small particles with negligible sedi- 
mentation velocity a TD = 1. 

In the present study we have assumed that the co- 
efficients Kd and Ktd are independent of the aerosol 
concentration because we consider the case m p N -C p, 
where m p is the particle mass. When this condition is not 
valid, nonlinear effects of aerosols on atmospheric turbu- 
lence should be taken into account. This is a subject of 
a separate study. 

Note that there is another effect that is called "tur- 
bophoresis" and related to the particle inertia. The tur- 
bophoresis results in an additional mean particle veloc- 
ity due to inhomogeneity of turbulence [see Caporaloni 
et al., 1975; Reeks, 1983]. The turbophoresis and tur- 
bulent thermal diffusion are totally different phenomena. 
Indeed, averaging Eq. (T2J) over fluctuations we obtain the 
mean particle velocity: 



(V p )< - U< - r, 



dXJj 
dt 



T s Vj{uiUj) + T s (uj(V r -u)) .(23) 



For example, in isotropic turbulence (uiUj) = 
(1/3) (u 2 )^-, and the mean particle velocity reads: 

V p = U- r s ^ ^V(u 2 ) +r s (u(V-u)) . (24) 
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The third term in Eq. (|24[) describes the effect of tur- 
bophoresis due to inhomogeneity of turbulence. The ra- 
tio of the mean particle velocity due to turbophoresis to 
the effective particle velocity caused by the phenomenon 
of turbulent thermal diffusion is 

| y t U rbo| _ ^ | V(u 2 ) | /(u 2 ) 

\Vf\ t T a TD \VT\/T ■ [ > 

Since t s -C t T , the mean particle velocity due to tur- 
bophoresis is much smaller than the effective particle ve- 
locity caused by the phenomenon of turbulent thermal 
diffusion. Although both effects are related to the parti- 
cle inertia, the effective particle velocity due to the phe- 
nomenon of turbulent thermal diffusion, originates from 
the turbulent particle flux (un), i.e., describes the col- 
lective statistical phenomenon, while the mean particle 
velocity due to turbophoresis originates directly from the 
expression for mean particle velocity. 

The mechanism of phenomenon of turbulent thermal 
diffusion is also principally different from molecular ther- 
mophoresis (and molecular thermal diffusion). The basic 
difference between these phenomena is explained in the 
following. The phenomenon of turbulent thermal diffu- 
sion occurs due to a combined action of turbulence effects 
and particle inertia effect, while molecular thermophore- 
sis is caused by purely kinetic effects related to thermal 
motion of molecules. In stratified turbulent flow turbu- 
lent thermal diffusion and molecular thermophoresis oc- 
curs simultaneously, although the effect of turbulent ther- 
mal diffusion for large Reynolds numbers is essentially 
stronger than the effect of molecular thermophoresis. In 
particular, the ratio of the effective velocity due to tur- 
bulent thermal diffusion to the velocity caused by molec- 
ular thermophoresis is of the order of Reynolds number. 
In this estimate we use formula for the molecular ther- 
mophoretic velocity V4h ~ ^|VT|/T [see Friedlander, 
2000] and Eq. (TT4"| for the effective velocity caused by 
turbulent thermal diffusion. For example, in the atmo- 
spheric turbulent boundary layer Reynolds numbers are 
of the order of 10 7 , that implies that molecular ther- 
mophoresis is negligibly small. 

Note that in the context of thermal convection or 
stably stratified turbulence, an anelastic approximation 
[V-(pu) = 0] is used for low Mach numbers. This is 
standard for describing deep convection but introduces 
the possibility of a compressible fluid velocity tied to 
variations in the fluid density. The fluid mass flux is 
divergence-free but a tracer particle moves with the fluid 
velocity. The flow compressibility is linked to the tem- 
perature fluctuations and there are dynamic correlations 
between fluid density, temperature and pressure. Irre- 
spective of particle inertia, this introduces a possibility 
of a mean drift even in homogeneous turbulence, i.e., 
causes turbulent thermal diffusion. 

In the next sections we discuss the potential impact of 
the phenomenon of turbulent thermal diffusion on dis- 
tribution of the atmospheric constituents. In particular, 
we analyze the observational information in the vicinity 



of strong temperature gradients (near the tropopause) in 
order to find out whether the effect of turbulent thermal 
diffusion is: (i) observable, (ii) significant, (iii) provides 
well-grounded explanations of the observed aerosol lay- 
ers and their positions. For the lower troposphere, we use 
the numerical simulations as a tool to study the contri- 
bution of the effect of turbulent thermal diffusion to the 
lower-troposphere vertical profiles of aerosol concentra- 
tion. Note also that the temperature minimum in the at- 
mosphere exists only for the absolute temperature while 
the potential temperature increases with altitude in the 
upper troposphere and the lower stratosphere. 



III. AEROSOL MEASUREMENTS BY GOMOS 
NEAR THE TROPOPAUSE 

A. Outline of the observational technique 

Tropopause is a well-known region with strong gradi- 
ents of temperature and also with substantial amount of 
particles, which remain there over long periods. Simulta- 
neous observations of the vertically-resolved aerosol con- 
centrations and temperature are scarce but the presence 
of the aerosol layers in the vicinity of the tropopause is 
well established [Brasseur and Solomon, 2005]. 

For joint analysis of temperature profiles and aerosol 
concentrations, we use a unique dataset obtained 
from GOMOS instrument (Global Ozone Monitor- 
ing by Occultation of Stars) onboard the Envisat 
satellite [Kyrola et al., 2004; Bertaux et al, 2004, 
http://envisat.esa. int/instruments/gomos|. GOMOS is 
equipped with the UV/Visible/NIR spectrometers, which 
record stellar spectra transmitted through the atmo- 
sphere continuously as the star sets behind the Earth 
limb. The measurements are performed in the limb- 
viewing geometry with the sampling frequency of 2 Hz. 
The atmospheric transmission spectra obtained after di- 
viding the stellar spectra observed through the Earth 
atmosphere by the reference spectrum, recorded above 
the atmosphere, contain spectral features of absorption 
and scattering by gases and particles. This allows recon- 
structing the vertical profiles of O3, NO2, NO3, O2, H2O 
and aerosol extinction in the atmosphere. The vertical 
sampling resolution of GOMOS data is 0.5-1.7 km. 

While ozone can be retrieved up to 100 km altitude, 
other species are usually detectable in the upper tropo- 
sphere and in the stratosphere. The lowest altitude of 
GOMOS measurements is from 5 km to 20 km; it de- 
pends mainly on stellar brightness and clouds top height. 
The stellar light can be transmitted only through thin 
clouds (like subvisual cirrus clouds), and they appear 
as increased aerosol extinction in GOMOS data. Usu- 
ally, aerosols are retrieved with a good accuracy down to 
tropopause and slightly below from occultations of bright 
stars. 

The GOMOS inversion of the chemical composition is 
performed in two steps [Kyrola et al., 1993]. First, atmo- 
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spheric transmission data from every tangent height are 
inverted to horizontal column densities (along the path 
of the light beam from the star) for gases and optical 
thickness for aerosols (spectral inversion). Then, for ev- 
ery constituent, the collection of the horizontal column 
densities at successive tangent heights is inverted to ver- 
tical density profiles (vertical inversion). 

Since the aerosol extinction spectrum is not known a 
priori, a second-degree polynomial model is used for the 
description of the aerosol extinction /3 aoro in the GOMOS 
retrievals: 

/3acro = y N(z) [1 + ci (A - A ref ) + c 2 (A - A rof ) 2 ] , 

(26) 

where A is wavelength in nm, A re f = 500 nm, a = 
3 x 10~ 7 cm 2 nm is the scaling factor, N(z) is the aerosol 
mean number density at altitude z (in the units of cm -3 ), 
and parameters c\ and c 2 determine the wavelength de- 
pendence of the aerosols extinction spectra. The aerosol 
number density N and the parameters c\ and c 2 are re- 
trieved from GOMOS data. 

High-resolution temperature profiles (HRTP) in the 
stratosphere and the upper troposphere are retrieved 
from the synchronous scintillation measurements by the 
GOMOS fast photometers operating at 1 kHz sampling 
frequency at red (650-700 nm) and blue (470-520 nm) 
wavelengths [Dalaudier et al., 2006]. The measurement 
principle exploits the chromatic refraction in the atmo- 
sphere. The bi-chromatic scintillations recorded by the 
photometers allow accurate determination of a refractive 
angle, which is proportional to the time delay between 
the photometer signals. The procedure of conversion of 
refractive angle profiles to temperature profiles is similar 
to that used in radio occultation. At altitudes 18-35 km, 
HRTP is retrieved with the vertical resolution 200-250 m 
and the accuracy of 1-2 K. The best accuracy is achieved 
in vertical (in orbital plane) occultations of bright stars 
[Sofieva et al., 2007]. Below 15 km, the quality of HRTP 
decreases due to low signal-to-noise ratio, broadening of 
scintillation peaks as a result of chromatic smoothing and 
the violation of the assumptions used in retrievals (in par- 
ticular, the weak scintillation assumption). 

B. Observed aerosol and temperature profiles 

For joint analysis of HRTP and aerosol profiles, two 
data-sets were selected: 217 successive occultations of 
one of the brightest stars, Canopus (referred hereafter 
as to S002), with visual magnitude -0.7 and the effec- 
tive temperature 7,000 K corresponding to location of 
the ray perigee point over 36 N-37 N (mid-latitudes) in 
Jan-Fcb 2003, and 560 successive occultations of S029 
(visual magnitude 1.6, the effective temperature 10,200 
K) located at 10-20 N, for the same period. The observa- 
tion period was selected arbitrarily. Random inspections 
of aerosol and temperature profiles at different locations 
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FIG. 2: Examples of aerosols (black lines) and temperature 
profiles (dashed grey lines) observed by GOMOS for mid- 
dle latitudes (36-37N) in January 2003, and the independent 
ECMWF temperature analysis data at measurement locations 
(grey solid lines) . 



and seasons at low and mid latitudes have indicated that 
the features described below are common. 

Since the vertical resolution of HRTP (250 m) is much 
finer than that of aerosol profiles (2 km), HRTP were 
smoothed down to resolution of aerosol profiles. For 
the sake of reliability, we also included the temperature 
profiles extracted from archive of operational analysis 
of European Centre of Medium-Range Weather Forecast 
(ECMWF). Usually, smoothed HRTP and ECMWF pro- 
files are in a very good agreement in the considered alti- 
tude range. 

It has been noticed that the aerosol concentration and 
temperature profiles are often anti-correlated (see a few 
examples in Figures 2 and 3). At a qualitative level, this 
is in agreement with the theoretical predictions related to 
the impact of the phenomenon of turbulent thermal dif- 
fusion on the aerosol profiles (the impact of other effects 
is considered in the discussion section). 

Similarities between the results in the low- and mid- 
latitude regions suggest a common mechanism behind the 
observed anti-correlation. The data for the equatorial 
tropopause, however, have to be treated with care: very 
low temperatures (T < 198 K) at the tropical tropopause 
are favorable for formation of cirrus clouds [Brasseur and 
Solomon, 2005], which are seen as increased aerosols ex- 
tinction in the GOMOS data. Nevertheless, current un- 
derstanding and observations of the cirrus clouds position 
them slightly below the tropopause, while in our cases the 
aerosol layer is practically symmetrical with regard to 
the temperature profile (unfortunately, the vertical reso- 
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S002, 36-37N, Jan-Feb 2003, 195 profiles S029, Jan-Feb 2003, 10-20N, 130 profiles 
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FIG. 4: Histogram of correlation coefficients between the 
ECMWF temperature profiles and GOMOS aerosol retrievals. 
Left: occultations of Canopus, mid latitudes (195 profiles out 
of 217 total covering the tropopause region), right: occulta- 
tions of S029, tropics (130 profiles out of 560 covering the 
tropopause region). For rejected profiles the retrievals have 
been terminated above the tropopause due to dense clouds. 




FIG. 3: As Figure 2, but for locations over tropics. All nota- 
tions are as in Figure 2. 



lution was insufficient for unambiguous conclusions). It 
must be noted that the cirrus cloud formation is not im- 
portant for mid-latitude profiles, because the observed 
temperature over 206 K is too high for the formation of 
the cirrus clouds. The latter usually occurs for the tem- 
perature smaller than 198 K (Jensen et al., 1996). 

The slope of aerosol extinction spectra can serve as 
an indicator of the aerosol type: large particle have flat 
spectra in UV/Visible range, while small particles have 
larger extinction in UV than in Visible range, which ap- 
proaches the A" 4 scattering law in case of very small par- 
ticles (see, e.g., discussion by Vanhellemont et al. [2005]). 
The analysis of the slope of the GOMOS aerosol extinc- 
tion spectra also indicates that the observed mid-latitude 
aerosols are probably background sulfate aerosols (small 
particles) rather than ice crystals (large particles). 

The statistically significant anti-correlation between 
the temperature and the aerosol concentration has been 
observed for 30 % profiles at mid-latitudes and for 50 % 
of profiles in tropics (see Figure 4). The distributions of 
the correlation coefficient between the relative gradients 
of temperature and the aerosol concentration are strongly 
skewed with modes close to —1. 

A further insight can be obtained from quantitative 
analysis of the profiles and, first of all, the functional 
relation between temperature and aerosol concentrations. 
The theory of turbulent thermal diffusion predicts the 
following relation between the steady-state profiles of the 
mean temperature and mean particle number density: 



N(z) [T(Z)} a TD = exp 



(27) 



which follows from Eq. (|2Tj) . Here the turbulent thermal 
diffusion ratio a TD is determined by Eq. (f2"2")) . Note that 
Eq. (|2T|) has been previously used by Eidelman et al. 
(2004, 2006a,b) for analysis of laboratory measurements 
of TTD. 

Equation (|27|) allows experimental determination of 
a TD from the measured profiles of temperature and 
aerosol concentration without using Eq. (|22j) . Indeed, 
assuming that the right-hand-side of Eq. (|27|) is constant 
above and below the temperature minimum (not neces- 
sarily the same constants) , yields the power-law type re- 
gression with unknown constant parameter. Taking loga- 
rithm and differentiating Eq. (|27| with respect to height, 
we arrive at a linear regression equation for determining 
a TD . Application of this procedure to GOMOS profiles 
with statistically significant negative correlation (at the 
significance level of 99 %) yields the histograms of a TD 
shown in Figure 5, which have a pronounced peak and a 
narrow width. As one can see in Figure 5, the values of 
a TD for 60-70 % of the anti-correlated profiles are very 
similar. This modal value, in turn, varies with latitudes 
being smaller for equatorial regions. 

In the comparison of the observed values of a TD with 
the theoretical predictions we have to take into account 
several parameters affecting it. The theoretical predic- 
tions for a TD [see Eq. (|2"2")) ] relate it to: 

(i) atmospheric pressure that affects kinematic viscos- 
ity and the mean free path of molecules; 

(ii) the coefficient of turbulent diffusion Kr>; 

(iii) aerosol size and material density. 

Theoretical dependencies of the turbulent thermal dif- 
fusion ratio a TD versus the coefficient of turbulent dif- 
fusion Kd are shown in Figs. 6-7, and the function a TD 
versus the size of aerosols is shown in Fig. 8. Theoret- 
ical values of a TD have been determined using Eq. (|22|) 
for Riy = and Lp = 8 km. Increasing Riy results in 
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S002, 36-37N, Jan-Feb 2003, 60 profiles S029, Jan-Feb 2003, 10-20N, 56 profiles 
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FIG. 5: Histogram of the TTD ratio a TD : a regression be- 
tween the relative gradients of temperature and number con- 
centrations of particles for profiles having statistically signifi- 
cant negative correlation. Left: mid latitudes (S002 measure- 
ments), right: tropics (S029 measurements). 
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6: Theoretical values of the TTD ratio a TD versus the 
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coefficient of turbulent diffusion Kd for the altitude 17 km 
and different sizes of aerosols: 5 pm (solid line), 2.5 pm 
(dashed-dotted line) and 1 pm (dashed line). Thin dashed- 
dotted line a TD — 1 corresponds to non-inertial particles. 



the increase of the turbulent thermal diffusion ratio a 



In calculations we used the following formulas for the 
Reynolds number, Re = 3Kr>/2C n i', and for the Stokes 
time t s = C c p a al / (18 p v), where C c is the Cunningham 
correction factor. 

Inspection of Figure 6 (the dependence of a TD at 17 
km altitude versus Kd and aerosol size) reveals that the 
variations of the turbulent thermal diffusion ratio a TD 
can be very large. The dashed line (aerosol of 1 pm 
in diameter) shows that a TD exceeds 7 when Kd ~ 1 
m 2 s _1 or less. The parameter a TD is much larger than 
those obtained in the laboratory studies [see Elperin et 
al., 2000b; Eidelman et al., 2004, 2006a, 2006b], which 
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FIG. 7: Theoretical values of the TTD ratio a TD versus the 
coefficient of turbulent diffusion Kd for aerosols of 1 pm size 
and different altitudes: 17 km (solid line) and 1 km (dashed- 
dotted line). Thin dashed-dotted line a TD = 1 corresponds 
to non-inertial particles. 
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18.1 km (solid line); Kd 
km (dotted line); K D = 3 x 10 5 cm 2 s 



to non-incrtial particles. 



have been conducted for much higher values of Kd and at 
surface pressure. The comparison with surface-pressure 
conditions for lpm aerosols (Figure 7) shows that lower 
air density at the tropopause almost doubles the parame- 
ter a TD , with another factor of 5-10 due to the lower eddy 
diffusivity. The combined impact of a lower fluid density 
and a smaller eddy diffusivity is illustrated in Figure 8, 
where the solid line corresponds to the tropopause con- 
ditions while other curves correspond to the case of the 
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lower troposphere. 

Comparison of the histograms in Figure 4 and Fig- 
ures 6-8 demonstrates that the observed values of a TD are 
fairly close to the predicted ones. The theory of turbu- 
lent thermal diffusion also provides a natural explanation 
for the lower values of a TD observed in tropics. Despite 
the higher altitude and the lower pressure, the dominat- 
ing parameter is the turbulent diffusion coefficient Kd, 
which is often high in the equatorial region [Parames- 
varan et al., 2003; Fujiwara et al., 2003; Yamamoto et 
al., 2003]. In fact, one can consider an inverse problem: 
the observed aerosol concentration and temperature pro- 
files can be used for determining the turbulent diffusion 
coefficient Kd in each particular case using a TD as an in- 
put parameter. This would yield the information about 
the level of turbulence near the tropopause obtained in- 
directly from the aerosol concentration and temperature 
profiles. 

In this section we have shown that the effect of tur- 
bulent thermal diffusion can be significant even for small 
particles. Indeed, for 1-2 /mi particles, the fluid velocity 
u is much larger than the incrtial term \T s d\i/dt\, and 
usually particles of this size are used in laboratory ex- 
periments as flow tracers. However, the phenomenon of 
turbulent thermal diffusion is determined by two con- 
tributions: (i) the correlation of the divergence of the 
inertial term, V- (r s du/dt), with turbulent velocity and 
(ii) the correlation of the divergence of the non-inertial 
term, V- u, with turbulent velocity. This results in effec- 
tive particle velocity that can be larger than the termi- 
nal fall velocity W g . Inspection of Figs. 6-7 shows that 
a TD can be much larger than 1 for Kd <C 10 6 cm 2 / s. 
Consequently, the contribution to the effective velocity 
of aerosols [see Eqs. (jT4"]) and (f2"2")) ] caused by the inertia 
term is larger than the contribution of the non-inertial 
term. On the other hand, for the lower troposphere (see 
the next Section), whereby the turbulent diffusion coeffi- 
cient Kd ~ (1 — 5) x 10 5 cm 2 / s, the parameter a TD w 1, 
and the effective velocity of aerosols is due to the non- 
inertial term. 



IV. ASSESSMENT OF THE CONTRIBUTION 
OF THE EFFECT OF TURBULENT THERMAL 
DIFFUSION TO THE LOWER TROPOSPHERE 
COMPOSITION 

Analysis of the role of turbulent thermal diffusion in 
the lower troposphere is quite involved due to highly dy- 
namic character of the atmospheric boundary layer, in 
particular, its strong diurnal cycle. Varying wind and nu- 
merous sources injecting aerosols at various heights pose 
additional difficulties. With these conditions, observa- 
tional evidence of the turbulent thermal diffusion in the 
lower troposphere is very difficult to detect, and a prelim- 
inary analysis of available radar data has not produced 
unequivocal results. 

The influence of turbulent thermal diffusion in the tro- 



posphere was evaluated via numerical modelling. This 
section discusses the means of including the phenomenon 
of turbulent thermal diffusion into the chemical transport 
models and presents results of several numerical exper- 
iments with the modelling system SILAM in order to 
reveal the impact of formulations of turbulent thermal 
diffusion on the predicted aerosol distribution. 



A. Setup of the modelling experiment 

The SILAM system [Sofiev et al, 2006] is a dual-core 
Lagrangian-Eulerian dispersion model. Its meteorologi- 
cal pre-processor provides the system with all parame- 
ters required for the simulations, including the input for 
TTD assessment. For the experiments, we used the Eulc- 
rian dynamic core based on original advection scheme by 
Galperin [2000] with the vertical turbulent diffusion pa- 
rameterized via K-theory and implemented in the model 
via extended resistive analogy by Sofiev [2002] . This ver- 
tical diffusion scheme meets the key requirement of the 
current experiment because it explicitly treats the parti- 
cle sedimentation and allows virtually any profile of ver- 
tical turbulent mixing coefficient Kd- 

The goal of the numerical experiment was to evalu- 
ate the impact of TTD on the aerosol concentrations 
in the lower troposphere. For that purpose, the model 
was run through a reference year 2000 over the Euro- 
pean continent using the archive of the limited-area me- 
teorological model HIRLAM [Unden et al., 2002]. As 
an aerosol tracer, we used primary anthropogenic par- 
ticles (i.e. aerosols directly emitted from anthropogenic 
sources) in two size classes: fine particles with less than 
2.5 /mi in diameter - PM 2.5, and a coarse-mode fraction 
from 2.5 /zm to 10 /zm size - PM 2.5-10. To simplify the 
simulations, both particle size ranges were considered as 
single bins with characteristic diameters of 1.4 /mi and 
6 /mi and densities of 1230 kg m~ 3 and 1500 kg m~ 3 , 
respectively. Exchange of particles between these bins 
due to condensation and coagulation was neglected. An- 
thropogenic emission for both classes was taken from the 
EMEP database (European Monitoring and Evaluation 
Programme, http://www.emep.int). 

The turbulent thermal diffusion velocity V° s was in- 
troduced as an additive term to particle sedimentation 
velocity W g so that the total mean vertical particle ve- 
locity reads: 



w g + vf = CcPaa * 9 



18/i 



K D ol„ 



V 2 T 
T 



(28) 



where z-axis is directed along the acceleration due to 
gravity g, a* is particle diameter, C c is the Cunningham 
correction factor and fi = pv is a dynamic viscosity of 
air. 

Vertical turbulent diffusion coefficient Kd is deter- 
mined using SILAM standard meteorological processing 
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routines, which evaluates Kd in a simplified form avoid- 
ing the uncertain parameters, such as turbulent kinetic 
energy or turbulent length scale. The computations of 
Kd start at the surface with evaluation of Kr>{z\ = 1 
m) after Gcnikhovich et al. [2004] and the height of the 
boundary layer Habl as described by Sofiev et al. [2006]. 
Following the surface-layer assumptions, it is assumed 
that Kd oc z in the range from z\ up to z\ = 0.1-ffABL, 
then Kd remains constant up to the height Habl- At 
the top of the boundary layer the coefficient K D sharply 
decreases by an order of magnitude and then remains 
constant up to the top of the modelled domain. These 
formulations are evidently too crude for detailed eval- 
uation of the near-tropopause phenomena discussed in 
the previous section but are sufficient for studying the 
lower troposphere processes, which were in the focus of 
the modelling experiment. 

The standard resistance-based diffusion term with tur- 
bulent and laminar-layer resistances represents the diffu- 
sion pathway of the aerosol dry deposition. For the stud- 
ied particle size range, in typical conditions and with- 
out turbulent thermal diffusion, the standard diffusion 
term is much smaller than the gravitational sedimenta- 
tion. However, in our simulations, this term has been 
included because the TTD velocity V° s can partly or en- 
tirely outweigh the sedimentation velocity, so that the 
diffusion component may become important. 



B. Results of the simulations 
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FIG. 9: Mean annual vertical profiles with (dashed lines) and 
without (black lines) TTD for PM 2.5-10 over two charac- 
teristic regions: Norwegian mountains (triangles) and Arctic 
Ocean (circles). Unit is fig m . 



FIG. 10: Mean annual near-surface concentrations of PM 2.5 
and PM 2.5-10 without the TTD term (panels a and b, re- 
spectively), and a ratio between the PM 2.5 and PM 2.5-10 
concentrations obtained with and without the TTD term in 
the SILAM model formulations (panels c and d, respectively). 
Unit of concentrations is /ig m -3 . 



The long-term averaging of the obtained results high- 
lighted the mean influence of the turbulent thermal dif- 
fusion: moderate in absolute values but systematic uplift 
of the aerosols. The driving force of this uplift is the 
decrease of temperature with increase of the altitude. 
Examples in Figure 9 for Norwegian mountains (panel 
a) and Arctic Ocean (panel b) show that, depending on 
the region elevation, the TTD impact can result in both 
increase and decrease of the near-surface concentrations. 
In the low-altitude regions and close to the sources, the 
near-surface concentrations decrease while over the ele- 
vated and mountain areas they grow (Figure 10). The 
pollution masses (aerosols), which are raised to higher 
altitudes by the TTD velocity V° s , appeared to be trans- 
ported over larger distances and at higher speed since the 
mean wind at the plume height is stronger and the dry 
deposition is weaker. Both effects increase the overall 
transport distance of aerosols. In turn, this resulted in 
wider distribution of aerosols, higher concentrations in 
remote areas and a larger fraction of mass leaving the 
source region, such as Europe. 

The above changes are moderate in absolute values. In 
our numerical experiment, the transport of coarse par- 
ticles PM 2.5-10 outside the modelled domain has in- 
creased by 5-15 % depending on season. For PM 2.5 
the difference was small. Episodically, the impact of 
TTD was more significant but the year-long simulations 



did not reveal any case when the daily concentrations 
changed by more than by the factor of 1.5 due to the 
TTD term. 

The difference between the particle size sections can 
be explained as follows (see Figures 6 and 7). In the 
lower atmosphere the turbulent diffusion coefficient is 
usually larger than 10 5 cm 2 s _1 . Therefore, for aerosols 
smaller than 2.5 /im the turbulent thermal diffusion ra- 
tio is a TD w 1 (as well as for gases), while for coarser 
aerosols the TTD ratio a TD varies from 5 to 10 (see Fig- 
ures 6 and 7). Therefore, the same temperature gradi- 
ent affects the dispersion of coarse aerosols 5-10 times 
stronger than that of the fine particles. This is also ev- 
ident from a simple fact that both velocities, and 
the sedimentation W g , are proportional to the squared 
particle size. 



V. DISCUSSION 

The phenomenon of turbulent thermal diffusion has 
been predicted theoretically for the atmospheric and lab- 
oratory turbulent flows and then detected in a series of 
laboratory experiments. Therefore, there is no doubt in 
the existence and importance of this effect in case of 
strong temperature gradients. The goal of this study 
is to assess the importance of this phenomenon in the 
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atmosphere and elucidate where one can find the atmo- 
spheric signature of this effect. The main problem is that, 
contrary to the controlled laboratory experiments, the in- 
fluence of TTD on aerosols distribution in the Earth at- 
mosphere has to be distinguished from other phenomena 
affecting aerosol profiles. 

The most important phenomena affecting the aerosol 
profile near the tropopause, are as follows: 

(i) aerosol dynamics; 

(ii) cloud microphysics; 

(iii) dynamic interaction of the vertical gradients of the 
turbulent diffusion coefficient Kd and aerosol concentra- 
tions; 

(iv) gravitational sedimentation of particles; 

(v) various dynamic phenomena violating the steady- 
state assumption and disturbing the aerosol profiles, e.g., 
tropical deep convection; 

(vi) variation of the aerosol concentration profiles 
caused by the Bernoulli effect and the wind speed gradi- 
ent near the tropopause. 

The contributions of these effects to the aerosol pro- 
file strongly vary and have to be evaluated case-by-case. 
Quantitative evaluation of impact of these effects often 
requires such details about aerosol composition and size 
spectrum, as well as about the atmospheric conditions, 
that are not available. However, some conclusions can be 
drawn from the presented data. 

There are several factors which indicate that the mech- 
anism related to TTD is universal and widespread, and 
should exist at least in low- and mid- latitudes. They 
include: 

(i) striking similarity of the aerosol profiles in tropics 
and mid-latitudes; 

(ii) very large fraction of considered aerosol profiles re- 
veal the temperature-concentration anti-correlation and 
faithfully follow the functional dependence given by 
Eq. ([27]). The theory of TTD provides the explanation 
for this anti-correlation and even predicts this functional 
dependence. 

In the subsequent qualitative analysis, we explore the 
possible alternative mechanisms which can maintain such 
aerosol profiles stable over a long time period. It must be 
emphasized that we do not consider the mechanisms of 
formation and transport of these amounts of particles to 
high altitudes. These can be associated with numerous 
effects, which are much faster and stronger then TTD. 
Our analysis is focussed on the mechanism which can ex- 
plain the existence the observed long-living aerosol con- 
centrations in different parts of the globe. 

Formation of the cirrus clouds can be of some impor- 
tance for tropical regions where the situation is quite 
ambiguous but it is hardly of importance for the mid- 
latitude GOMOS profiles. Notably, the aerosol optical 
characteristics are closer to sulphates than to water ice 
crystals [see Vanhellemont et al., 2005]. Therefore, cirrus 
clouds formation cannot be the main mechanism respon- 
sible for the observed aerosol profiles. It should be also 
kept in mind that the GOMOS observations are feasible 



only in the absence of thick clouds which, consequently, 
cannot affect the aerosol concentration profiles. 

The dynamic interaction of aerosol concentration 
near the tropopause with the fast-changing Kd with 
altitude follows from expression: V z (-?TdV z N) = 
(V ' Z K D )(V Z N) + K D V 2 Z N. The first term in the r.h.s. 
of this equation corresponds to advection of the species 
in the direction opposite to the gradient of the diffusion 
coefficient. This would lead to temporal accumulation 
of the particles near the altitude with sharp changes of 
the turbulent diffusion coefficient Kd, which is closely 
related to local temperature gradients. However, this ex- 
planation also has several limitations. 

Firstly, the temperature minimum at the tropopause 
exists only for the absolute temperature while the poten- 
tial temperature is steadily growing with altitude in the 
upper troposphere and the lower stratosphere. Conse- 
quently, the stratification above and below the minimum 
of the temperature is stable, which implies Kd is com- 
paratively small in the considered height range and its 
gradient is also not large. 

Secondly, Eqs. ([5]) and © in their "classical" form with 
V z s = do not have a steady-state solution with a max- 
imum of concentration inside a modelled domain. The 
increased concentration at the top of turbulent domain 
can occur only in case of strong elevated aerosol sources 
and will not remain for a long time. Eventually, this 
maximum will be smeared by turbulence and finally will 
disappear. Sedimentation velocity, that is substantial 
at these altitudes due to low air density and viscosity, 
would speed-up this process. However, numerous ran- 
domly chosen GOMOS observations (in addition to the 
ones reported above) show that the aerosol layers in the 
vicinity of the temperature minimum are commonly oc- 
curring. 

In general, the strong correlation between the temper- 
ature and aerosol concentration profiles cannot be ex- 
plained by any of the dynamic processes whose main 
driving forces are associated with the wind and turbu- 
lence rather than temperature per-se. Therefore, such 
mechanisms cannot cause the observed functional depen- 
dence between temperature and aerosol concentrations. 
Hence, the effects, such as Bernoulli jet, are also hardly 
of primary importance here. 

Formation of new aerosols and deep convection in trop- 
ical regions are probably the most-prominent of the re- 
maining alternative mechanisms of the aerosol appear- 
ance at high altitudes. However, neither of them ex- 
plains the aerosol persistence and formation of the pro- 
files anti-correlated with temperature. Indeed, sedimen- 
tation of particles is substantial at these altitudes, and 
such downward flux of particles would be readily observ- 
able and would result in essentially different forms of the 
relation (|2"7|) for the parts of the aerosol profiles above 
and below the temperature minimum. Since nothing of 
this kind is observed, we have to assume that some other 
mechanism is responsible for formation of the vertical 
particle profile after aerosols have been formed or trans- 
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ported by, e.g., deep convection towards the tropopause. 

The above qualitative analysis shows that none of the 
considered effects can explain the observed stable in time 
widespread functional relation between the temperature 
and aerosol concentrations given by Eq. (|27p with a TD 
given by Eq. (|22|) (see numerous examples in Figures 6- 
8). Validity of these relations for the analyzed hun- 
dreds of cases is probably the main argument in favour of 
the phenomenon of turbulent thermal diffusion as one of 
the main mechanisms that stabilizes the observed anti- 
correlation between the temperature and aerosol con- 
centration profiles, and renders this anti-correlation a 
widespread feature of aerosol profiles at low and middle 
latitudes. 

It is worth mentioning that the near-tropopause condi- 
tions are drastically different from the laboratory exper- 
iments [see Eidelman et al, 2004; 2006a, 2006b; Buch- 
holtz et al., 2004]. At normal pressure, the typical val- 
ues of a TD measured in the laboratory experiments do 
not exceed 3-4, while the near-tropopause observations 
showed much higher values. The main reason for this 
difference is substantially lower air density and viscosity, 
as well as the smaller turbulence intensity. With corre- 
sponding corrections taken into account, the theoretical 
predictions for a TD appear to be in a fairly good agree- 
ment with the observations (see Figures 6-8). 

Detection of the turbulent thermal diffusion in the 
observations in the atmospheric boundary layer is even 
more complicated than in the middle atmosphere. Our 
modelling experiments showed that characteristic times 
of the competing effects in the boundary layer are signif- 
icantly shorter than those required for approaching the 
steady-state solution with the evident signature of the 
turbulent thermal diffusion. In addition, the TTD ratio 
GLr-p jrj IS substantially lower near the surface than in the 
upper troposphere, thus requiring stronger temperature 
gradients for the effect to be easily detectable. Conse- 
quently, a comprehensive analysis of profile-resolving ob- 
servations is required for a detection and quantification 
of the phenomenon of turbulent thermal diffusion in the 
lower troposphere. 

The numerical experiments with the SILAM model 
provided the first estimates for the TTD impact in 
the lower troposphere and highlighted its main features. 
Moderate but systematic changes in the aerosol distribu- 
tions that emerged from the simulations seem to comply 
well with the theory of this phenomenon. Indeed, TTD 
affects the whole mass of aerosols usually acting oppo- 
site to gravitational sedimentation and being comparable 
with it. For some naturally occurring temperature gra- 
dients and intensity of turbulence, sedimentation of fine 
particles can be fully compensated by V£ , essentially 
canceling their dry deposition (the diffusion-driven dry 
deposition is usually insignificant for coarse particles). In 
our numerical simulations, the average effect of the tur- 
bulent thermal diffusion was at the level of 10-15 % with 
general tendency of redistribution of aerosol masses up- 
wards. As an immediate consequence, the aerosol trans- 



port distance has increased as well as the concentrations 
over the elevated and remote regions (see Figure 10). 

As can be seen from the maps shown in Figure 10, 
the relative effect of TTD on coarse particles is consid- 
erably stronger than that on fine aerosols. Indeed, the 
TTD term [see Eqs. @, {H} and [fig])] has a contribu- 
tion which is proportional to the squared particle size. 
Consequently, for particles with small sedimentation ve- 
locity W g the absolute impact of TTD is also small and 
easily overshadowed by other effects, such as macro-scale 
turbulence, transport with mean vertical wind, etc. For 
coarser particles, the sedimentation is much more impor- 
tant, therefore its adjustment due to TTD is also more 
pronounced. 

Notably, Eq. (QJJ) for the effective velocity, V° s , for 
for non-inertial particles or gaseous admixtures (i.e., for 
a TD = 1) has been already applied by Atreya et al. [1999] 
for study of vertical mixing in the atmospheres of Jupiter 
and Saturn. This equation was obtained by Atreya et al. 
[1999] using phenomenological arguments. 

The theory of the phenomenon of turbulent thermal 
diffusion still needs to be refined by more tests in carefully 
controlled experiments, atmospheric observations and in 
direct numerical simulations of density-driven turbulent 
flow, where the flow physics can be controlled and mod- 
ified to test the different aspects of the theory. 



VI. CONCLUSIONS 

This study is a follow-up of theoretical investigations 
and laboratory experiments [Elperin et al., 1996; 1997a; 
Eidelman et al., 2004, 2006a, 2006b], which demonstrated 
the existence of phenomenon of turbulent thermal diffu- 
sion. This paper provides the first observational evidence 
and quantitative evaluation of the importance of TTD in 
the atmosphere. 

The analytical part of this paper provides a set of prac- 
tically applicable formulations, directly useable for anal- 
ysis of observations and ready for implementation in dis- 
persion models. The most important effect is the extra 
term in the diffusion equation expressed through macro- 
scale parameters of turbulence and aerosol features, see 
Eqs. ©, GD-(IS) and o. 

Application of these formulations to the analysis of 
GOMOS observations near the tropopause seems to ex- 
plain otherwise confusing shape of aerosol vertical profiles 
with elevated concentrations located almost symmetri- 
cally to the temperature profile. To the best knowledge 
of the authors, such symmetry has not been explained 
before. 

It must be emphasized that the characteristic times 
and the magnitude of turbulent thermal diffusion are 
insufficient for formation of aerosol layers. Numerical 
simulations have confirmed that the impact of TTD is 
systematic but moderate, which also complicates detec- 
tion of this phenomenon among dynamic processes in 
the boundary layer. Hence, turbulent thermal diffusion 
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alone is insufficient for elevating substantial amount of 
aerosols to the tropopause and formation of the aerosol 
layers with the enhanced concentration of particles. In 
the performed modelling experiment, the simulation time 
and/or the size of the computational domain were insuf- 
ficient for the formation of aerosol concentrations. Con- 
sequently, it is plausible to expect that turbulent thermal 
diffusion hardly plays the main role in the formation of 
the near-tropopause aerosol layers. Other effects, such 
as deep convection, buoyant plumes from wild-land fires, 
formation of new aerosols, etc., may be more important 
in this respect. However, turbulent thermal diffusion is, 
probably, the only effect that can trap the aerosols in the 
vicinity of the minimum of temperature once they are 
formed. 



APPENDIX A: THE EFFECTIVE VELOCITY 
FOR SMALL PECLET NUMBERS 

In this Appendix we derive Eqs. ([6]) and (fT0| using 
a quasi-linear approach or a second order correlation ap- 
proximation [see, e.g., Moffatt, 1978]. We rewrite Eq. ([5]) 
in the following form: 



^ + V-Q-DV 2 n 

at 



I 



(Al) 



where Q = vn — (vn) is the nonlinear term and / = 
AV-v + (v-V)A is the source term. Let us neglect 
the nonlinear term but keep the molecular diffusion term 
in Eq. (|Alj) . For this reason this approach is called a 
quasi-linear or perturbational approach. This approxi- 
mation for a given velocity field is valid only for small 
Peclet numbers (Pe <C 1), where Pe = uq£/D. Let us 
rewrite Eq. (|A1|) in a Fourier space. Then the solution 
of Eq. (|Aip is given by n(uj, k) = Gd(uj, k)/(oj, k), where 
G D (u,k) = {Dk 2 + ito)- 1 . 

We apply a standard two-scale approach, i.e., the non- 
instantaneous two-point second-order correlation func- 
tion can be written as follows: 



(Ui(ti,x)n(i 2 ,y)) = / («i(wi,k 1 )n(a;2,k2))exp[i(k 1 -x 
+k 2 -y) + i(w±ti + W2*a)] du>i duj 2 dki dk 2 
= J F^ n) (uj, k) cxp[ik-r + iu f] duo dk , (A2) 

where we use large scale variables: R = (x + y)/2, K = 
ki + k 2 , t = (ti + t 2 )/2, = UJ-y + u>2, and small scale 
variables: r = x — y, k = (ki — k2)/2, f = ti — t 2 , 
uj = (ui\ — uj 2 )/2 [see, e.g., Roberts and Soward, 1975]. 
Here 

(uj, k) = J (vi(uji, ki) n(oj 2 , k 2 )) exp[iflt 

+iK-n]dndK, (A3) 

and w\ = uj + £1/2, ui 2 = -uj + £1/2, ky = k + K/2, 
k 2 = — k + K/2. We assume here that there exists a 



separation of scales, i.e., the maximum scale of random 
motions I is much smaller than the characteristic scales 
of inhomogeneities of the mean particle number density 
and mean fluid density. 

Hereafter we use the simplest model for the second 
moment of a random velocity field in a Fourier space: 



(ui(uj, k) Uj(— uj, — k)) 



<E{k) 



8irk 2 



5{uj) , (A4) 



that satisfies the continuity equation in anelastic approx- 
imation V u = UjAi, where Aj = —Vip/p, S(uj) is the 
Dirac's delta-function, Sij is the Kronecker tensor, the 
energy spectrum function is E(k) = k^ 1 (q— 1) {k/k a )~ q , 
the exponent 1 < q < 3, the wave number fc = l/l, the 
length £ is the maximum scale of random motions and 
Mo is the characteristic velocity in the maximum scale of 
random motions. Hereafter we neglect the small terms 

We consider non-inertial particles, i.e., v — u. Then 
the turbulent flux of particles is given by 

F[ n) = J (ui(uj,k)I(-uj,-k)}G D (-uj,-k)du>dk.(A5) 

After integration in u-space and in k-space we arrive at 
Eq. ^ for the turbulent flux of particles with 



(q-l)u £ 

Kd = ^nr Pe ' 

Vf = K D ^P = -K D ViT 



(A6) 
(A7) 



P T 
In the derivation we used the integral: 

sin Odd (ki k 3 /k 2 ) dip = (4tt/3) ^ 

(i 



APPENDIX B: THE EFFECTIVE VELOCITY 
FOR LARGE PECLET NUMBERS 

In this Appendix we derive Eqs. ((6]) and (|T0|) using the 
t approach that is valid for large Peclet and Reynolds 
numbers. Using Eq. (|A1|) written in a Fourier space we 
derive equation for the instantaneous two-point second- 
order correlation function (t, k) = (ui(t, k) n(t, — k)): 



dF 



{n) 



= ( Ui (t,k)I(t,-k)) +MF { l III \k) , (Bl) 



where MF % {III) (k) = -[(u.V-Q) + ((d Ui /dt)n) - 
D(ui V 2 n)]k are the third-order moment terms appear- 
ing due to the nonlinear terms which include also molec- 
ular diffusion term. 

The equation for the second moment includes the first- 
order spatial differential operators M. applied to the 
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third-order moments F^ Iir K A problem arises how to 
close the system, i.e., how to express the third-order 
terms AiF^ 111 ^ through the lower moments F^ 11 ^ [see, 
e.g., Orszag, 1970; Monin and Yaglom, 1975; McComb 
1990]. We use the spectral r approximation which pos- 
tulates that the deviations of the third- moment terms, 
M.F( II]C '(\s), from the contributions to these terms af- 
forded by the background turbulence, MF^ 111 ' ^^), can 
be expressed through the similar deviations of the second 
moments, 

MF (III \k) - MF<- ni >°\k) 

1 "Hk) 



.(fc) 



(B2) 



[see, e.g., Orszag, 1970; Pouquet et al., 1976; Elperin et 
al., 2006], where r r (k) is the scale-dependent relaxation 
time, which can be identified with the correlation time 
r(fc) of the turbulent velocity field for large Reynolds 
and Peclet numbers. The functions with the superscript 
(0) correspond to the background turbulence with a zero 
gradient of fluid density and a zero gradient of the mean 
number density of particles. Validation of the r approxi- 
mation for different situations has been performed in nu- 
merous numerical simulations and analytical studies [see, 
e.g., review by Brandenburg and Subramanian, 2005; and 
also discussion in Rogachevskii and Kleeorin, 2007, Sec. 
6]. 

Note that the contributions of the terms with the su- 
perscript (0) vanish because when the gradients of fluid 
density and the gradients of the mean number density 
are zero, the flux of particles vanishes. This implies that 



Eq. jB2| reduces to MF^'^k) = -F} n> (k)/r(k). We 
also assume that the characteristic time of variation of 
the second moment -F^(k) is substantially larger than 
the correlation time r(k) for all turbulence scales. There- 
fore, in a steady-state Eq. (|B1[) yields the turbulent flux 
of particles F^ = J r(k){Ui(t, k) I(t, -k)) dk. Now we 
use the following simple model for the second moment of 
turbulent velocity field: 



(tt,(k) Uj(-k)) 



\E{k) 



8wk 2 



-- (Aikj - A,/ 



After integration in k-space we arrive at Eq. 
turbulent flux of particles with 

P 



K off = K 



D " 



= -Kr 



T 



(B3) 
for the 

(B4) 
(B5) 



In the derivation we used the following expression for the 
turbulent correlation time, r(fc) = 2tr (k/ko) 1 ~ q , where 
ty = £/uq is the characteristic turbulent time. Com- 
parison of Eqs. (|A7j> and (|B5|) shows that the form of 
the effective velocity T^ off is the same for small and large 
Peclet numbers, while the values of turbulent diffusion 
coefficients are different in these two cases. 



APPENDIX C: THE BUDGET EQUATION FOR 
THE TURBULENT FLUX OF PARTICLES 



Let us consider inertial particles suspended in the tur- 
bulent fluid. Their concentration is characterized by the 
mean value, N, and fluctuation, n, of the number density 
of particles [m -3 ]. The budget equation for the turbulent 

flux of particles, F^ — (vt n), reads 



DF 



in) 



Dt 



-(u 4 V-v)N - e 



3 
in) 



(CI) 



where D /Dt = d/dt + U-V, v is the particle velocity, 
u and U are the fluctuations and mean components of 

(n) 

the fluid velocity, <frj ■ is the third-order turbulent flux 



of F- 



e\ is the molecular dissipation rate of F^ : 



in). 



in) 



(ui w l n) + - (p n) 6ij , 



.(«) 



-D (ui V 2 n) -v(n V 2 u t ) , 



(C2) 



(C3) 



v is the kinematic viscosity of fluid, D is the Brownian 
diffusivity of particles and p are the fluid pressure fluctu- 
ations. The Kolmogorov closure hypothesis implies that 



.in) 



F 



in) 



C n tT 



pin) E l/2 



(C4) 



where l z is the vertical turbulent length scale, E z is the 

1/2 

vertical turbulent kinetic energy, = £ Z /E Z is the tur- 
bulent time, C' n is an empirical dimensionlcss coefficient. 
The term in Eq. (JCI) is given by the following ex- 
pression: 



)(n) 



1 



{pVi-n) + (3ei (n0 p ) , 



(C5) 



p is the potential temperature fluctuation, e is the ver- 
tical unit vector, (3 = g/T* is the buoyancy parameter, g 
is the acceleration of gravity, and is a reference value 
of the mean temperature T and p is the fluid density. In 
Appendix D we derive the following expression for the 



(«). 



correlation term Q 



Q\ n) = ~C D e t 



+F ( " ) v,e), 



E 



1/2 



{F j V j N+{6 p V-v)N 



(C6) 



where Cd is an empirical dimensionless constant, Fi — 
(ui 6p) is the flux of the potential temperature and 
is the mean potential temperature (or the mean virtual 
potential temperature accounting for specific humidity), 
that is defined as 9 = T(P,/P) 1 -^" 1 . Here P is the 
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mean pressure, P* is its reference value, and 7 = c p /c v = 
1.41 is the specific heat ratio. Then Eq. (|C1[) can be 
rewritten 



DF 



00 



Dt 



Vi$i n) = -r ij V j N-Fj n) V j U i 



F {n) E l/2 



2El' 2 



C n l z 

x [(e i F j +e j F i )\7 j N + 2e l (9 p V-v)N 
+ (e^+e^V.e] , (C7) 

where the r.h.s. describes the local budget of the tur- 
bulent flux of particles (aerosols), F^ n \ The third-order 

(n) 

flux <f>z2 can be expressed using the turbulent diffusion 
approximation: 



= -C F KV Z F^ , 



(C8) 



1/2 

where K = £ z E z . Then the budget equation for the 
vertical particle flux, F z n \ reads: 



DF, 



00 



Dt 



C F W Z (KV Z F^ 



N 



T s (U Z V 2 p) 



f3F z e z 



Et 



- (**■-<*'-&)( 



(C D (3 



V Z N — N 



y z p 



v 2 e 



E r 



(;„) 



CnW El 12 



(C9) 



where t s is the particle Stokes time and Cp is the em- 
pirical dimcnsionlcss constant. Equation (|C9[) is a com- 
plementary equation to the non-local closure model sug- 
gested by Zilitinkevich et al. [2007]. In the steady-state, 
homogeneous regime of turbulence, Eq. (|C9|) yields the 
vertical component of the turbulent flux of particles in 
the form of Eq. ([6]) with the coefficients Kd and Ktd 
given by Eqs. (fT5]) and p^|) . Equation (|C7|) also allows 
us to determine the horizontal components of the particle 
flux [see Eq. (20)]. 

Now let us derive Eq. (TB"]) . The effective velocity V 
depends on the turbulent heat flux (u6) that is deter- 
mined by the well known equation: (u6) = — VT 
[see, e.g., Monin and Yaglom, 1975], where Kh ~ 
is the coefficient of turbulent thermal conductivity. Note 
that herein we do not consider situation with very high 
gradients when gradient transport assumption is vio- 
lated. The above formula for the mean turbulent heat 
flux is written in the r-space. The corresponding sec- 
ond moment in k-space is given by (ltj(k) 8{— k)) = 
—r(k) (ttj(k) Uj(— k))Vj-T, where the second moment of 
turbulent velocity field is given by Eq. (|B3|) . the turbu- 
lent correlation time is r(k) — 2 tr (fc/fco) -2 ^ 3 , the energy 
spectrum function is E(k) = (2/3)fc _1 (fc/fc ) -5 / 3 . 

Multiplying equation for (ui(k) 6(— k)) by 
—k 2 r(k) and integrating in k-space we obtain 
i T (uV 2 6>) = (2/3) ln(Re) J B(Re,a») VT, where the 
function i?(Re, a*) = 1 when the particle size a* < a cr , 



and _B(Re,et„) = 1 — 3 ln(a*/a cr )/ ln(Re) when a* > a cr 
[see Elperin et al. , 2000a] . Here a cr is the critical particle 

size that is given by a cx = £ v (p/ ' Pp) 1 ^ 2 an d t-u = ^Re~ 3 ^ 4 
is the Kolmogorov viscous scale of turbulence. 



APPENDIX D: DERIVATION OF THE BUDGET 
EQUATION FOR {n6 p ) AND EQ. (fC6]) 

Equations for fluctuations of the potential temperature 
and the particle number density read: 



-UjVj (Q + 9 p ) + kV 2 6 p , 



(Dl) 



DOp 
Dt 
Dn 

— = -{ujVj + V-v) (JV + n) - n V- V + D V 2 n , 

(D2) 

where V is the mean particle velocity. Multiplying 
Eq. (|D1|) by n and Eq. (1D2|) by p , averaging and adding 
yield the budget equation for the correlation function 

E„a = (nd n ): 



V-v)N-E n0 V-V-e (ne) , (D3) 



Here is the third-order turbulent flux: 

$W=<«i»0p>> (D4) 
and e( n8 " > is the molecular dissipation rate of E n g: 



T (n9) 



-D(9 p V 2 n) ~ K (nV 2 9 p ) 



(D5) 



Molecular dissipation rate can be expressed using 

the Kolmogorov closure hypothesis: 



(n6) En8 E n g E, 



1/2 



C n gtT C n g l z 



(D6) 



where C n g is the dimcnsionlcss constant. 

In the steady-state, homogeneous regime of turbulence, 
Eq. (|D3p reduces to the turbulent diffusion formulation: 



E nB = -Cng—^iF^VjG + FjVjN 



e z 

9 P V-v)iV) , 



(D7) 



where we consider only gradient approximation neglect- 
ing higher spatial derivatives. 

(n) 

Now let us determine the correlation term Q\ 
which includes the correlation function p _1 (pVjn) (see 
Eq. l|C5p). To this end we use the following identities: 

p _1 {fl P v»p) = -i3{e p A- x wle p ) , (D8) 
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which follow from the Navier-Stokes equation. Indeed, 
calculating the divergence of the Navier-Stokes yields: 



v 2 P = -pv 



(D9) 



Applying the inverse Laplacian to Eq. (|D9|) we arrive at 
Eq. (IDS)) . Therefore, 



P -1 (pV<n) = ^((Vi^A-Vi^) =/3V i (nA- 1 V^ p ) 
-^(«A^V z V,fl p ), (D10) 



Our analysis shows that the last term in Eq. (|D10[) can 
be estimated as: 



1 - 



arctan y/5 



(Dll) 



[see Zilitenkevich et al., 2007], where S = and l z and 
lh are the correlation lengths of the correlation function 
(n(t, x) 9 p (t, y)) in the vertical and horizontal directions. 
For nearly isotropic case (8 <C 1), Eq. (|D1 1|1 reads: 



,A" X V? 
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En 



(D12) 



while for a strongly anisotropic case (6 ^> 1), Eq. (|D1 1|) 
yields: 

(n A" 1 V^ p ) « (l - JL + . (D13) 



Equations (|D7|1 - (|D13|) yield the expression for the cor- 
relation term Q^ n ' = (pVjTi) + /3ej (n0 p ) given by 
Eq. (IC61) . 
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